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1. Introduction 

Since the discovery of D 4 (2317) [jl]] by BABAR in April 2003, a series of new heavy mesons 
with open-charm and closed-charm have been observed by Belle, CDF, CLEO, BABAR, and BES. 
Among these new heavy mesons, the most intriguing ones are the charmonium-like states, X (3872) 

7(3940) 01, 7(4260) 0], Z(3930) [§], and X(3940) [@]. Evidently, one can hardly interpret 
all of them as orbital and/or radial excitations in the charmonium spectrum. Thus it is likely that 
some of them are exotic (non-qq) mesons (e.g., molecule, diquark-antidiquark, and hybrid meson). 
Theoretically, the central question is whether the spectrum of QCD possesses these resonances, 
with the correct quantum numbers, masses, and decay widths. Now the most viable approach to 
tackle this problem from the first principles of QCD is the lattice QCD with exact chiral symmetry 

©SB H. 

Recently, we have investigated the mass spectrum of closed-charm exotic mesons with J PC = 
1 [11] and 1 ++ j ji2[ ] respectively, in the framework of lattice QCD with exact chiral symme- 



try. By constructing molecular and diquark-antidiquark operators with quark content (cqcq), we 
measure the time-correlation function of each operator, and fit to the usual formula 

W { e -rnat + e -ma{T-t)} 

to extract the mass (m) and the spectral weight (W) of the lowest-lying state overlapping with the 
exotic meson operator. From the ratio of the spectral weights of two spatial volumes (20 3 and 
24 3 ) with the same lattice spacing, we can infer whether this state is a resonance (W20/W24 ~ 1) or 
2-particle scattering state (W20/W24 ~ (24/20) 3 = 1.728). 

Our results suggest that 7(4260) and X(3872) are in the spectrum of QCD, with J PC = 1 
and 1 ++ respectively, and both with quark content (cucu). Note that we have been working in the 



isospin limit (m u = m^), thus our results [11, 12] cannot exclude the possibility of the existence of 



exotic mesons with quark content (cdcd), even though we cannot determine their mass differences 
from those with (cucu). Moreover, we also observe heavier exotic mesons with quark contents 



(cscs) and (cccc), for J PC = 1 ++ [|12[|, and 1 [|lTj] respectively. 



Now if the spectrum of QCD does possess exotic mesons with quark content (cqcq), then 
it is likely that there are also exotic mesons with other quark contents, e.g., (cqcs) and (cscq). 
However, in general, whether any combination of two quarks and two antiquarks can emerge as 
a hadronic state relies on the nonperturbative dynamics between these four quarks. Recently, we 
have investigated the lowest-lying mass spectrum of molecular and diquark-antidiquark operators 



with quark content (cscq)/(cqcs) and J = 1 + [13]. Our results suggest that there exists a J = 1 



resonance around 4010 db 50 MeV. It is interesting to see whether this state will be observed by 
high energy experiments. 

2. Lattice setup 

To implement exact chiral symmetry on the lattice, we consider the optimal domain-wall 



fermion proposed by Chiu [14]. From the generating functional for «-point Green's function of 



the quark fields, the valence quark propagator in background gauge field can be derived as [ 14 ] 



(q(x)q(y)) = (D c + m q ) x ^, D c = 2m j * ^J^l i 2A } 
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where m q is the bare quark mass, mo is a parameter in the range (0,2), S(H W ) is equal to the 
Zolotarev approximation of the sign function of H w (H w = JsD w , and D w is the standard Wilson 
Dirac operator minus mo). In the limit N s — ► °° (where N s + 2 is the number of sites in the 5- 
th dimension), D c is exactly chirally symmetric, i.e., D c 75 + YsD c = 0. In the continuum limit, 
the valence quark propagator (D c + m q )~ l — > [/^(^ + z'gA^) + m^]~ 1 . A salient feature of the 
valence quark propagator ( |2~I| ) is that the bare quark mass m q is well-defined for any gauge field 
configuration, unlike the Wilson quark propagator. In practice, we compute the valence quark 
propagator by the nested conjugate gradient 

D{m q )Z = {m q + (mo-m q /2)[l + Y 5 S{H w ))}Z = l (2.2) 

S(H W )Z = Y^bizu (Hl + ci)zi = Z (2.3) 
l 



where b\ and q are Zolotarev coefficients [15|, |16Q . To attain the maximal efficiency of our system 



(CPU time and memory usage vs. the precision of the sign function S(H W )), we use Neuberger's 



2-pass algorithm [|17| |18p in the inner conjugate gradient (|2.3|). 

We generate 100 gauge configurations with single plaquette gauge action at j3 = 6.1, for two 
lattice volumes 24 3 x 48 and 20 3 x 40, with the same lattice spacing. Fixing mo = 1.3, we project 
out 16 low-lying eigenmodes of \H W \ and perform the nested conjugate gradient in the complement 
of the vector space spanned by these eigenmodes. For N s = 128, the weights {co s } are fixed with 
Knin = 0.18 and X max = 6.3, where < X(\H W \) < Xm ax for all gauge configurations. For each 
configuration, point-to-point quark propagators are computed for 30 bare quark masses in the range 
0.03 < m q a < 0.8, with stopping criteria 10~ n and 2 x 10 -12 for the outer and inner conjugate 
gradient loops respectively. Then the norm of the residual vector of each column of the quark 
propagator is less than 2 x 10 



-11 



\\[p c + m q )Y- I|| < 2 x 10" 11 , Y = {1 -m q /{2m )}- l {Z- l/(2m )}, 
and the chiral symmetry breaking due to finite Af s is less than 10 , 

Z f S 2 Z 



ztz 



1 



< 10 



-14 



for every iteration of the nested conjugate gradient. 

Then we measure the time-correlation functions of pseudoscalar and vector meson operators. 
We determine the inverse lattice spacing aT l = 2.237(76) GeV from the pion time-correlation 
function, with the experimental input of pion decay constant f n = 13 1 MeV. The strange quark bare 
mass m s a = 0.08 and the charm quark bare mass m c a = 0.80 are fixed such that the corresponding 
masses extracted from the vector meson correlation function agree with 0(1020) and //i/a(3097) 
respectively [ 191 2C ] . 



3. Four-quark meson operators 

In QCD, there is no reason why any hadronic state composed of two quarks and two antiquarks 
cannot emerge as a resonance. Thus, it is interesting to see whether such four-quark mesons can 
exist in the spectrum of QCD, and also to identify them with any meson states observed in high 
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energy experiments. To tackle this problem in the framework of lattice QCD, one needs to construct 
interpolating operators which can have good overlap with the 4-quark meson states. Otherwise, 
their signals may not be unambiguously identified, due to the limitations in statistics, the lattice 
size, as well as the unphysically-heavy u(d) quark masses which have to be chirally extrapolated 
to their physical values. Further, for any lattice calculations in quenched approximation, one should 
avoid any operators (e.g., scalar operator cq) which have strong quenched artifacts as m q — > m u ^. 

In general, there are two ways to construct 4-quark meson operators: (i) molecular operator, 
and (ii) diquark-antidiquark operator. In the following, we only present the 4-quark operators which 
have good overlap (except F4) with the hadronic states in our investigations [JTTJ, [T^J, [13|]. Explicitly, 
they are 



T PC 



V 



jPC = 1+H 



where 



M 3 (x) = -^={(qy 5 yc) x (cy 5 q) x - (cy 5 tfqMqy 5 c).*} (3-1) 

Y 4 (x) = -^{(q T Cy 5 yc) ax (qCy 5 c T ) ax + (qCy 5 y[c T ) ax (q T Cy 5 c) ax } (3.2) 

Mi(x) = -^{(qyc) x (cy 5 q) x - (c^q) z (qy 5 c) x } (3.3) 

X 4 (x) = -±= {(q T Cyc) ax (qCy 5 c T ) ax - {qCyJc T ) ax {q T Cy 5 c) ax ) (3.4) 

M s (x) = -j= {(qyc) x (cy 5 s) x - (cys) x (qy 5 c)^} (3.5) 

D 4 (x) = { (q T Cric) ax (sCr 5 c T )a, - (sCylc T ) ax (q T Cy 5 c\ LX } (3.6) 



(qrQ).* — qaaxT a pQ a p x 



denotes a meson operator with antiquark field q coupling to quark field Q through the Dirac matrix 
r. Here x, {a,b,c} and {ot,j8} denote the lattice site, color, and Dirac indices respectively. The 
diquark operator is denoted by 

(q T TQ) ax = SabcqbaxTapQcpx 

where e a i, c is the completely antisymmetric tensor, C is the charge conjugation operator satisfying 
Cy^C -1 = — y^ and (Cys) r = — Cy^. Thus the diquark transforms like color anti-triplet. For T = 
C/5, it transforms like J p = + , while for T = Cy (i = 1,2,3), it transforms like 1 + . 



4. Results and discussions 

In Table [T[ we summarize our results of the masses of the lowest-lying states of the 4-quark 
operators obtained in Refs. [|ll], 12, 13], where in each case, the first error is statistical, and the 
second one is our estimate of combined systematic uncertainty including: (i) possible plateaus (fit 
ranges) with ^ 2 /d.o.f. < 1; (ii) the uncertainties in the strange quark mass and the charm quark 
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mass; (iii) chiral extrapolation (for the entries containing u/d quarks); and (iv) finite size effects (by 
comparing results of two lattice sizes). Note that we cannot estimate the discretization error since 
we have been working with one lattice spacing. 

4.1 J PC = 1 states 

The molecular operator M3 ~ {(q75%c)(c75q) — (c757q)(q75c)} detects a resonance with mass 
4238(31)(57) MeV in the limit m q -> m H>d , which is naturally identifed with 7(4260). This sug- 
gests that 7(4260) is indeed in the spectrum of QCD, with quark content (cucu) and J PC = 1 . 

For the diquark-antidiquark operator Y4 ~ {[q r C75yc] [qC7sC r ] + [q r C75c] [qCysjfc 7 ]}, it also 
detects a state with mass 4267 (68) (83) MeV, in the limit m q — ► m u . We suspect that this state might 
be the same resonance captured by the molecular operator M3. However, we are not sure that this 
state is a resonance since the ratio of spectral weights (R = W20/W24) for two different lattice 
volumes with the same lattice spacing deviates from one (the criterion for a resonance) with large 
errors as m q — > m u . It is plausible that such a deviation is due to the quenched artifacts which can 
be evaded if one incorporates internal quark loops. 

Now, in the quenched approximation, our results suggest that 7(4260) has a better overlap 
with the molecular operator M3 than the diquark-antidiquark operator Y4. Whether this implies that 
7(4260) behaves more likely as a D\D molecule than a diquark-antidiquark meson is subjected to 
further investigations, especially those incorporating dynamical quarks. 

For molecular and diquark-antidiquark operators with quark fields (cscs), they both detect a 
resonance around 4450 ± 100 MeV, and for the molecular operator with (cccc), it detects a reso- 
nance around 6400 ± 50 MeV. These serve as predictions from lattice QCD. 

4.2 J PC = 1++ states 

Both the molecular operator M) and the diquark-antidiquark operator X4 detect a 1 ++ reso- 
nance around 3890 ± 30 MeV in the limit m q — > m u ^, which is naturally identified with Z(3872). 
This shows that X(3872) is indeed in the spectrum of QCD, with quark content (cucu), and 
J PC = 1 ++ . It is interesting to see that Z(3872) has good overlap with the molecular operator 
Mi as well as the diquark-antidiquark operator X4. This is in contrast to the case of 7(4260), in 
which 7(4260) seems to have better overlap with the molecular operator M3 than any diquark- 
antidiquark operators. This seems to suggest that X(3872) is more tightly bounded than 7(4260). 
It would be interesting to see whether this picture persists even for unquenched QCD. 

For m q = m s , heavier exotic meson resonance with J PC = 1 ++ is also detected, with quark 
content (cscs) around 4100 ±50 MeV. This serves as a prediction from lattice QCD. 

4.3 J PC = 1+ states 

Both the molecular operator M s and the diquark-antidiquark operator D4 detect a 1 + resonance 
around 4010 ± 50 MeV in the limit m q — > m^. Since its mass is just slightly above Y(3940), high 
energy experiments should be able to see whether such a resonance, say, X s , exists in some decay 
channels, e.g., X s — ► KnJ/^i 1 , in the near future. 
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J PC (J P ) 


Operator 


Mass (MeV) 


R/S 


Candidate 




-^lYuy^y-cVcysu) — Ccy<;'K-u')fuy<;c')l 


4238(31)(57) 


R 


Y(4260) 


I"" 


-4=[fsysy;c)(cyss) — (cysys)(sysc)l 


4405(3 1)(44) 


R 




j — 


(cY;c)(cc) 


6411(25)(43) 


R 




1 — 


^ { [u r Cy 5 yc] [uCy 5 c T ] + [u r Cy 5 c] [uCysfc 7 ] } 


4267(68)(83) 


R? 


Y(4260) 


1 — 


j- { [s r Cy 5 7ic] [sCy 5 c T ] + [s T Cy 5 c] [sCy 5 y ( r c r ] } 


4449(40)(55) 


R 




1 + + 


4=[(uyc)(cy 5 u) - (cy ; -u)(uy 5 c)] 


3895(27)(35) 


R 


X(3872) 


1 + + 


^[(syc)(cy5s)-(cys)(sy 5 c)] 


4109(21)(32) 


R 




1 + + 


^ {(u r Cyc)(QCy 5 c T ) - (u r Cy 5 c)(uCyf c 7 )} 


3891(17)(21) 


R 


X(3872) 


1 + + 


^ {(s r Cyc)(sCy 5 c r ) - (s r Cy 5 c)(sCyf c r )} 


4134(19)(25) 


R 




| + 


^[(uyc)(cy 5 s) - (cys)(uy 5 c)] 


4007(34)(31) 


R 






j- {(u r Cyc)(sCy 5 c T ) - (u r Cy 5 c)(sCy r c r )} 


4015(25)(27) 


R 





Table 1: Masses of the lowest-lying states of the molecular operators and the diquark-antidiquark operators. 
The column R/S denotes resonance (R) or scattering (S) states. The entries with empty slots in the last 
column are theoretical predictions [O, O, O ], 



5. Concluding remarks 



Note that we are working in the quenched approximation which in principle is unphysical. 
However, our previous results on charmed baryon masses [19], charmed meson masses, and our 



predictions of /o and /d [ ]2Q| ] are all in good agreement with the experimental values [ 21| , [22| , |23 



24]. This seems to suggest that it is plausible to use the quenched lattice QCD with exact chiral 



symmetry to investigate the mass spectra of the charmoniumdike 4-quark meson operators, as a 
first step toward the unquenched calculations. The systematic error due to quenching can only be 
determined after we can repeat the same calculation with unquenched gauge configurations. (Note 
that Monte Carlo simulations of unquenched gauge configurations in lattice QCD with exact chiral 
symmetry, on the 20 3 x 40 and 24 3 x 48 lattices at j3 = 6.1, still remains a challenge to the lattice 
community.) However, we suspect that the quenched approximation only results a few percent 
systematic error in the mass spectrum, especially for charmed mesons. Apparently, to gauge how 
well the quenched approximation in determining the mass spectrum of 4-quark meson operators is 
to see whether the predicted states in Table |l| (entries with empty slots in the last column) will be 
observed in high energy experiments. If it turns out that high energy experiments do not find any of 
these predicted states, then it might imply that: (i) the mass spectra of 4-quark mesons in quenched 
QCD could be dramatically different from those in unquenched QCD, or (ii) the constancy of the 
spectral weight for different volumes with the same lattice spacing (i.e., WvJWv 2 — 1) might not 
be sufficient to guarantee that it is a resonance, or both (i) and (ii). On the other hand, if the 
predicted states are confirmed by high energy experiments (even if our predicted masses turn out 
to be 5 ~ 8% off the experimental values), it would be a remarkable success for (quenched) lattice 
QCD with exact chiral symmetry. Obviously, no matter what is the outcome, it will be interesting. 
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